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Abstract. Back-to-Back Correlations (BBC) of particle-antiparticle pairs are predicted to appear 
■ if hot and dense hadronic matter is formed in high energy nucleus-nucleus collisions. The BBC 

are related to in-medium mass-modification and squeezing of the quanta involved. Although the 
suppression of finite emission times were already known, the effects of finite system sizes and of 
collective phenomena had not been studied yet. Thus, for testing the survival and magnitude of the 
effect in more realistic situations, we study the BBC when mass-modification occurs in a finite sized, 
thermalized medium, considering a non-relativistically expanding fireball with finite emission time, 
QO ' an d evaluating the width of the back-to-back correlation function. We show that the BBC signal 

\Q , indeed survives the expansion and flow effects, with sufficient magnitude to be observed at RHIC. 
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In the late 1990's the Back-to-Back Correlations (BBC) between boson-antiboson 
pairs were shown to exist if the particles masses were modified in a hot and dense 
mediumyj]. Not much longer after that, it was also shown that a similar BBC existed 
between fermion-antifermion pairs with medium-modified masses [2]. A similar formal- 
ism is applicable to both BBC cases, related to the Bogoliubov-Valatin transformations 
of in-medium and asymptotic operators. Both the bosonic (bBBC) and the fermionic 
(fBBC) Back-to-Back Correlations are positive and have unlimited magnitude, thus dif- 
fering from the identical-particle correlations, also known as HBT (Hanbury Brown & 
Twiss) correlations, which are limited for both cases, being negative in the fermionic 
sector. BBC were expected to be significant for pj < 2 GeV/c. Nevertheless, already 
in the Ref.Ll], it has been shown that the duration of the emission process significantly 
suppresses its magnitude. The effects of finite system sizes and of collective phenomena 
were not known, which motivated us to investigate their consequences. Some prelimi- 
nary results will be discussed here and illustrated in some particular cases. 

0.1. Infinite Homogeneous Medium 

In our analysis, we assume the validity of local thermalization and hydrodynamics, 
from the beginning up to the system freeze-out, as well as a short duration of the particle 
emission. We also consider that the effective in-medium Hamiltonian can be written as 



H = H - / dxJy<Kx)SM 2 (x-y)<Ky), 



(1) 



where 

H = Jdx(<p 2 + \V<t)\ 2 + m 2 (l> 2 ) (2) 

is the asymptotic (free) Hamiltonian in the matter rest frame, and the second term 
in Eq. (UJ) describes the medium modifications. The scalar field represents quasi- 
particles propagating with momentum-dependent medium-modified mass m*, related to 
the vacuum mass, m, by 

m 2 (|k|)=m 2 -5M 2 (|k|). 

This implies that the dispersion relations in the vacuum and in-medium are given, 
respectively, by C0% = m 2 + k 2 and Q 2 = m 2 + k 2 = CO 2 - <5M 2 (|k|), where Q. is the 
frequency of the in-medium mode with momentum k. 

The annihilation (creation) operator, b (b^), for the in-medium, thermalized quasi- 
particles is related to the annihilation (creation) operator of the asymptotic, observed 
quanta with momentum W 1 = (fflfc,k), a (a'), by the Bogoliubov-Valatin transformation 



c k b l + s -k b -k ; cik = c k b k + s*_ k tf_ kl (3) 



where c k = cosh(/fc) and s k = sinh(/^); f k = jlog(^) is the so-called squeezing param- 
eter, since the Bogoliubov transformation is equivalent to a squeezing operation. 
On the other hand, the two-particle probability distribution is given by 



alal 2 a k2 a kl ) = (a J k a kl ) (4 2 «k 2 ) + (ala kl )(a k a kl ) + {a J k a J k2 ) (a kl a kl ) 



(4) 



The first term on the r.h.s. is proportional to the product of the two single-inclusive dis- 
tributions, with momenta i.e., N\(ki) = (0k,^ = ft>k, (^k^k,-), while the second term 

is proportional to the square modulus of G c (ki,k2) = G c (l,2) = V /C0 kl o\ 2 (a ki a kl ), 
the so-called chaotic amplitude. The last term is related to this new contribution, 
Gi(ki,k2) = Gj(l,2) = ycOk, o\ 2 {a k{ a kl ), which is called squeezed amplitude. 

In cases where the particle is its own anti-particle (for or 00 boson pairs, for 
instance), both terms contribute, and the full correlation function is written as 

r(k k 1 1 i l G C ( k i> k2 )l 2 i |G,(ki,k 2 )| 2 

C 2l K!,K 2 J f Gc(ki>kl)Gc(k2)k2) + Gc(kl)kl ) Gc(k2!k2 )' W 

where the first two terms correspond to the HBT correlation, and last term, represents 
this additional contribution to the correlation function, i.e., the squeezing part. 

In the above equations, the amplitudes were written in terms of the asymptotic oper- 
ators. However, the averages in Eq. © are estimated by means of the density operator 
p, as the thermal average for globally thermalized gas of b quanta that is homogeneous 
in the system with a certain volume V, with p = (1/Z) exp[— (V/r)(l/(27r) 3 ) / QJrb]. 
Being so, the above equations should be expressed in terms of the in-medium operators 
by means of Eq. ©, prior to performing the thermal averages. 

We know that the maximum value of the HBT correlation function is attained when 
ki = &2 = k, resulting in C c (k,k) = 2. Accordingly, the maximum value of the BBC 



correlation function is attained for k\ = —ki = k, resulting, after performing the thermal 
averages, into 

Qk,-k =1+ k 7?f — , (6) 

m(k)m(-k) 

where «i(k) = [|ck| 2 nk + |*-k| 2 («-k + 1)] is related to the spectral function by iVi (k) = 
(0^n\ (k); «k is the Bose-Einstein distribution function of the in-medium quanta 

with energy D. k at temperature T. Strictly speaking Eq. d6|) is valid only in the rest frame 
of the medium. 



0.2. Finite size medium moving with collective velocity 

For a hydrodynamical ensemble, both the chaotic and the squeezed amplitudes, 
G c and G s , respectively, can be written in the special form derived by Makhlin and 
Sinyukov [3] (see Eqs. (22) and (23) of Ref. [1]), namely 



G c (ki,k 2 



d 4 o^ 
d 4 o^(x 



I 1 Jiix x 



■K\ 2 e 



ki,2| 2 «i,2 (x) + \s-i - 2 \ 2 [n_i _ 2 (x) + I 



(7) 



G,(k l7 k 2 )= J -^W^ 2e 2^ 1>2 .x| / _ i;2C2 _ in _ i;2W+Cij _ 2/ _ 2;i [ni _ 2(x) + 1] |. 



(8) 



In Eq. © and ©, do^(x) = J 3 E j[i (x, z)F(x)dz is the product of the normal-oriented 
volume element depending parametrically on x (freeze-out hyper- surface parameter) and 
on its invariant distribution, F(x); a' 1 is the hydrodynamical freeze-out surface. 

For studying the expansion of the system we adopt the non-relativistic hydrodynam- 
ical model of Ref. ||4|]. In this model the fireball expands in a spherically symmetric 
manner with a local flow vector given by the four- velocity = y(l,v), assumed to be 



non-relativistic, with y = (1 — v 



- 1 / 2 « l+v 2 /2, where 
v=(u)r/R, 



being (u) and i? the mean expansion velocity and the radius of the fireball, respectively. 
We then divide the inhomogeneous medium into independent cells and assume that the 
expressions for G c and G s can be evaluated locally within each cell using Eq. ©. The 
squeezing coefficient can be written, in more detail, as 
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log 



K*](x)u v (x) 
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log 



(o ki (x) + co kj (x) 
£l ki (x)+Cl kj {x) 



f±U±M) 



(9) 



where, as in HBT, the pair momentum difference and the pair average momentum 



In 



are given, respectively, by q^j(x) = kf(x) —k^(x), and K^j(x) = ■. 

addition, we consider the Boltzmann limit of the Bose-Einstein distribution for n^, i.e. 
rii j ~ exp [— {Kf jU^ — n(x))/T(x)], and assume a time-dependent parametric solution 

of the hydrodynamical equations, i.e., fi{x)/T{x) = Hq/Tq — r 2 /(2R 2 ), as in Ref. 0]. 





FIGURE 1. Schematic representation of the region where the mass-shift occurs: on the l.h.s., the 
modification is extended to the whole system, whereas on the r.h.s. it happens only in a smaller portion of 
the thermalized medium. The figures represent cross-sectional areas of the Gaussian profiles. 



Furthermore, we consider two possible scenarios for the freeze-out: the first, corresponds 
to a sudden freeze-out, in which F(x) 5(t — To). The second scenario corresponds to a 
smeared freeze-out, for which ^L_B>2 ^-^-^M^ This last more realistic scenario has a 
dramatic effect on the Back-to-Back Correlation function, as already showed in Ref.Jl]], 
by reducing severely the signal's magnitude, even for a smearing of about At ~ 2 fm/c. 

In discussing finite-size effects, we distinguish between the volume of the entire 
thermalized medium, denoted by V (with radius R), and the volume filled with mass- 
shifted quanta, denoted by V s (with radius R s ). Naturally, V s < V in the general case. In 
the derivation of the expressions for G c ( 1 , 2) and G s ( 1 , 2) , for simplicity, we introduce a 
Gaussian profile function in the integrands, i.e., ~ exp [— r/ (2R) ]. In Fig. 1 we illustrate 
this by showing cross-sectional areas corresponding to Gaussian profiles, for the cases 
with V = V s and V > V s . 

In the non-relativistic limit, the accounting for the squeezing effects can be simplified 
for small mass shifts (m* — m)/m<^ 1, such that the squeezing parameter in Eq. © can 



be written simply as f(i, j, r) « ^ log [j^j ■ This limit is important, because in this case 

the coordinate dependence enters the squeezing parameter / only through the possible 
position dependence of the mass-shift which, in principle, could be calculated from 
thermal field models in the local density approximation. Therefore, in an approximation 
such that the position dependence of the in-medium mass can be neglected, the c(i, j) = 
cq and s(i, j) = sq factors can be removed from the integrands in Eq. © and © and all 
what remains to be done are Fourier transforms of Gaussian functions. 

For completeness, we write below the expression of the Back-to-Back Correlation 
function for the case where the mass shift occurs in entire volume of the system, V. 
A detailed discussion and more complete formulation of the problem, including the 
case of mass-shift in a smaller portion of the system, can be found in Ref.o]. In 
what follows, we will concentrate on the value of momenta of the participant pair that 
maximizes the BBC signal, i.e., the case in which ki = — k2 = k, using the fact that 
the single-inclusive distribution depends only on the absolute value of the momentum, 
i.e., G c (k, k) = G c (— k, — k). Nevertheless, since the strict condition of back-to-back 
pair holds only in the rest frame of the medium, we implicitly are considering the case 
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FIGURE 2. In the left panel, parts (a) and (b) illustrate the effect of finite emission times. The dashed 
curves, corresponding to a sudden emission (At = 0), have been decreased by a factor of 400, and the solid 
curves show the suppression caused by a finite emission duration, of about At ~ 2 fm/c, which drastically 
reduces the BBC magnitude. In the right panel we illustrate the cases with and without flow, for two 
values of the momentum of the back-to-back pair (in (c), k = 500 MeV/c and in (d), k = 1000 MeV/c). 
The mass-shifting is supposed to occur in the entire system volume. In the plots, m* is the in-medium 
modified mass and T stands for the freeze-out temperature. 



corresponding to a weak flow coupling in the expanding system, which is expected to be 
fulfilled for pair emitted near the center of the system. In this case, the BBC correlation 
function can then be written as 



C BBc(k, -k) = l + |c s 

R 3 {\c \ 2 + \so\ 2 ) 



2R 3 



l + : 
k 2 



m 2 <u> 2 ^- 3/2 



exp 



T 



2m+T 



1 + 



K" m 2 <u > 2 k 2 

ni2 u : ^ 3/2 ex P [ - y " 2^f + 2m 2 T 2 (\ + ^4^) 



+# 3 M 2 



(10) 



In Figure 2 we illustrate some of the results found in this non-relativistic approxima- 
tion, in the particular case of weak flow coupling. We see that the cases corresponding 
to the absence of flow and to its inclusion produce similar results within the limits of 
our illustration. However, depending on the value of k\ = —k2 = k, there are noticeable 
differences. Being so, we see that, for smaller values of k, the presence of flow seems to 
slightly enhance the signal, whereas at large values of k, the non-flow case wins. Nev- 
ertheless, the non-flow case grows faster with increasing k. In this particular example 
discussed here, the effect is mainly due to the denominator in Eq. (flOb . since it contains 
the flow parameter in the exponential's coefficient. 



0.3. Summary and Conclusions 



Our main goal on presenting these new results here was to revive the discussion on the 
search of the squeezed BBC. For fulfilling this purpose, we estimated the strength of the 
BBC signal in a more realistic situation, considering the mass-shifting in a finite region, 
and the emission occurring during a finite time interval. We also considered that the 
system expands non-relativistically and analyzed the simplified situation of weak flow 
dependence (i.e., back-to-back pairs emitted close to central system region) of the BBC. 
Finally, we employed 00 pair correlation for illustration. We showed in Figure 2 the 
back-to-back correlation function versus the in-medium shifted mass, m*. We also saw 
that both the non-flow and the flow cases produced similar results, with pronounced 
maxima around m fa m*. Although we did not show here results corresponding to 
the case of mass-shift occurring only in a smaller portion of the system, it is shown 
elsewhere||5l|g] to decrease proportionally to the size of the mass-shift region. However, 
the effect of decreasing the system size is far less significant than the finite emission 
time for reducing the BBC magnitude. We also saw that, at least in its weak-coupling 
limit, the flow may work for slightly enhance the BBC signal for small values of 
the momentum k. Our main conclusion, however, is that, in the particular framework 
discussed above, a sizeable strength of the squeezed BBC signal could be seen, making 
it a promising effect to be searched for experimentally at RHIC. 

Naturally, for a more refined calculation, it is mandatory to introduce a model-based 
mass-shift. After that, it is also essential to perform more realistic calculations with flow, 
in a less particular kinematic region, while simultaneously searching for those which 
could optimize the observation of the BBC signal. Also an estimate of the shape and 
width of the BBC around the direction k\ = —k2 = k should be implemented. Finally, for 
being able to make predictions closer to the experimental conditions, it will be extremely 
important to obtain some feed-back on the experimental acceptance, conditions, and 
restrictions that could finally lead to the BBC discovery. 
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